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We show, using the quench action approach T|, that the whole post-quench time evolution of 
an integrable system in the thermodynamic limit can be computed with a minimal set of data 
which are encoded in what we denote the generalized single-particle overlap coefficient Sjf°(A). This 
function can be extracted from the thermodynamically leading part of the overlaps between the 
eigenstates of the model and the initial state. For a generic global quench the shape of s*°(A) in 
the low momentum limit directly gives the exponent for the power law decay to the effective steady 
state. As an example we compute the time evolution of the static density-density correlation in 
the interacting Lieb-Liniger gas after a quench from a Bose-Einstein condensate. This shows an 
approach to equilibrium with power law t~ 3 which turns out to be independent of the post-quench 
interaction and of the considered observable. 

PACS numbers: 02.30.Ik,05.70.Ln,75.10.Jm 


Introduction. Understanding the non-equilibrium 
time evolution of a many-body interacting system is 
one of the main challenges in contemporary physics 
Hug. The study of systems with nontrivial interactions 
among their constituents is hard enough when the 
system is in its ground state; things however get even 
more complicated out of equilibrium, since most of the 
usual theoretical tools then become inapplicable. This is 
mainly due to the high energy regions of the spectrum 
which are probed by the time evolution, where the 
mean field approach and low energy approximations are 
not valid. Numerical simulations on the other hand 
are severely limited in the range of time or of system 
sizes g]. Developing new methods able to predict the 
short-, intermediate- and long-time dynamics of an 
out-of-equilibrium system is thus an urgent priority, 
especially in view of the rapid progress achieved in 
experiments on ultracold atoms [SHE]. 

Since the beginnings of quantum mechanics m , much 
interest has been devoted to the fundamental problem of 
calculating the time dependence of physical observables 
in states which are not eigenstates of the Hamiltonian 
driving the time evolution. This situation has now come 
to be known as a quantum quench mm and has been 
of major interest both from experimental and theoretical 
points of view. Most of the theoretical research focused 
so far on the expectation values of local observables at 
late times after the quench, when the system is in an ef¬ 
fective steady state. In particular the Generalized Gibbs 
Ensemble (GGE) hypothesis [TS1HB] focuses on the possi¬ 
bility of reducing the huge complexity of the initial wave 
function to a reduced set of information, incorporated in 
the local conserved quantities of the system, which gives 
the expectation values of all physical observables in the 
steady state. However the question of how to perform 
an analogous simplification for the whole post-quench 


time evolution, much more relevant from the experimen¬ 
tal point of view, has been poorly addressed, except in a 
few cases m nn iijjszi - 

The quench action method introduced in pP has 
recently proved to offer a procedure whereby one can de¬ 
rive, from first principles, not only the steady state itself 
but also the actual time evolution of physical observables 
Q 221 ES I58H35] . In summary, this method combines 
knowledge of initial state overlaps with functional 
integration techniques to extract the thermodynamically 
relevant information on the relaxation dynamics of an 
integrable system. The purpose of this Letter is to 
show that this approach, combined with the recent 
observations on the structure of the overlaps between 
eigenstates of different Hamiltonians [38[ 144H47] , is able 
to provide the full post-quench time evolution in terms 
of a reduced set of data which can be extracted from 
the thermodynamically leading part of the overlaps. 
It turns out that the same function, the generalized 
single-particle overlap coefficient Sg°(A), fixes the 
steady-state expectation values and the whole time 
evolution from t = 0 + after the quench. This is treated 
analogously to a system at thermal equilibrium with a 
sub-entropic gas of independent particle-hole excitations 
around the steady state constituting the whole effective 
spectrum necessary to compute the time dependence of 
all physical observables. A restricted class of excitations 
is then clearly seen to be the most relevant for the 
long-time behavior, giving a picture reminiscent of a 
field theory description of the asymptotics of correlations 
in equilibrium situations |48H50j . 

This letter is organized as follows. First we show how, 
for a generic integrable model, the quench action method 
[1] allows to extract the whole post-quench time evolu¬ 
tion from the complex function s^°(A) denoted here as 
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the generalized single-particle overlap coefficient. Then 
we specialize to the time evolution of the static density 
moment 32(3-’ = 0 ,t) of the interacting Lieb-Liniger gas 
after a quench from the ground state of the bosonic free 
theory. The same quench has been studied in a number 
of recent works HSIESIEDEBI. We use here the exact 
results for the post-quench saddle point reported in f5B] 
and we obtain a rare full post-quench time evolution of 
a physical observable in a truly interacting model that is 
closely related to recent experiments 1531531. 

Time evolution in an integrable model We consider 
an initial state I’I'o) which is not an eigenstate of the 
one-dimensional integrable Hamiltonian H for N parti¬ 
cles moving on a system size L with periodic boundary 
conditions. In a generic integrable model each eigen¬ 
state is specified by a set of N quantum numbers I = 
{Ij}jLi- The set of nonlinear coupled Bethe equations 
maps the mutually excluding quantum numbers in N 
quasi-momenta, called rapidites, I —► A = {Ay } jW x which 
take value in the complex plane. These are related to the 
one-particle momentum fco(A) and the scattering phase of 
the model 6{ A) [591 

2ttI N 

k 0 (X t ) = £ 0(A 4 - A fc ) i = l,...,N. (1) 

k =1 

All the possible different choices of quantum numbers 
I give a complete basis of eigenstates |A) with energy 
IS [A] which allows, given a local operator O, to resolve 
the time evolution of its expectation value on the initial 
state ('F o\e iHt de - iHt \^ 0 ) = (6(t)> 

A /i. 

( 2 ) 

where we introduced the overlap coefficients S 1 ^ 0 between 

the initial state and the eigenstates (H/ 0 1 A) = e _s A°. The 
double sum in ([2]) can be performed in general when the 
number N of constituents of the system is small. How¬ 
ever one is in general interested in the thermodynamic 
limit limth = limjv,L-xx> with fixed density n = N/L. 
The quench action approach introduced in [T] allows to 
move from a sum over the discrete representation for the 
eigenstates, in terms of the quantum numbers I to a func¬ 
tional integral over smooth distributions of rapidites and 
simple excitations over them. Given a smooth function 
p(X) of rapidities on the real axis with its normalization 
given by the density of particles f_ dX p(X) = n (under 
the string hypothesis it can be generalized to complex 
rapidities [60]), there is an entropic number ~ e s vv[p] 0 f 
finite size states that share the same expectation values 
of local operators [59]. The entropy is given by m 



( 3 ) 


where the density of holes is given by the total density 
p h = p* — related to the density of particles by the 
Bethe equations Q in the thermodynamic limit 

P* = ^ (fco + 6 * p) , (4) 

where we introduce the convolution between two func¬ 
tions f * g = dpf(X — p)g(g) and the derivative 
respect to A, = /. For later convenience we also 
introduce the scalar product on the real axis / • g = 
dp.f(p,)g(p). After restricting to the appropriate 
sub-Hilbert space with nonzero overlap (if discrete sym¬ 
metries are present), the overlaps become a smooth func¬ 
tional over the eigenstates. In particular they can be 
written as an extensive universal part (dependent only 
on the distribution p) with subleading corrections which 
depend on the finite number of particle-hole excitations 
{hi,Pi }™-li over the distribution p (which corresponds to 
displacing a number m of quantum numbers of one of the 
finite size state |A) — > \p ) which discretizes the distribu¬ 
tion p) 

St° -+ S[p) + 5s[p , {h^PijZx} + 0(1/N) , (5) 

where both quantities are given in terms of the general¬ 
ized one-particle overlap coefficient Sq°(A) 

S[p) =Ls*°-p, 

m 

5s[p, {huPi}^} = (so°(Pk) - So°(h k ) - F k ■ s* 0 ) . 

k =1 

( 6 ) 


The back-flow F k ( A) for a single particle-hole is computed 
in terms of the distribution p 

2irF k ^ =d(X~p k )~d(X-h k ) + e*F k . (7) 
P 


Therefore in the thermodynamic limit, for any weak op¬ 
erator (62j . we can write its time-dependent expectation 
value © as 


hm{6(t)) = 

th 



— 25R«S [p]+Sy-r [p] 


00 p 

Y / d [h,P)m 

m =0 ^ 


e -Ss[ P Ah,,P,}T=i]-iM P Ah,,P,}T=d^ p \d\ Pi {h l:Pl }^ 


+ mirr , 


( 8 ) 


with / d[h,p] m = UT=i f-o o dhjPihj) JY dpjp h (pj) 
denoting the sum over the macroscopic particle-hole 
excitations and mirr indicating the same sum as in ([8]) 
but with excitations on the left state [63]. The energy 
of a state 2£[A] —*■ £\p] + oj[p, {hi,Pi} r ^L^\ is given in 
terms of the one-particle energy eo(A) analogously to the 
overlaps © 

£[p\ = Le 0 - P, 

m 

uj[p, {hi,Pi}™, 1 ] = Y ( £ o(Pfc) - £0 (hk) - F k ■ e 0 ) ■ 
k =1 
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The matrix elements (p\0\p) can be computed by choos¬ 
ing one of the possible (large) finite size realizations 


|A) —► |p) of the distribution p( A) and using = 

(■*■ + • The same can be done for the off 

diagonal ones. Given these ingredients the sum in ([2]) 
can be evaluated in the saddle point ss s ^ = 0 of 

P P=Psp 


the quench action S®[p] = 25ft<S[p] — Syy[p] leading to 
an expression for the whole post-quench time evolution 
in the thermodynamic limit [TJ I38j 


1 . , I" 
lim(0(t)) = 2 

Z m=0 ^ 

0 -Ss[p sp ,{hi, Pi }l 




x e 


{Psp\0\psp,{hi,pi}™ =1 ) 


+ mirr. 

(9) 


A notable consequence of formula is that all the infor¬ 
mation to reconstruct the entire post-quench time evo¬ 
lution is contained in the function s^°(A) which can be 
extracted by taking the scaling limit of the overlap coef¬ 
ficients 


N 

St° = 55(4 0 (A J ) + O(^v- 1 )) . (10) 

3 =1 

The behavior of the exponent of s^°(A) around A ~ 0 
determines the power law for the large time relaxation of 
any physical observable. In the limit of large t we can 
indeed approximate the sum in ([9]) with the contribution 
of the saddle point itself and of the single particle-hole 
excitations 

1 f 00 

lim(0(t)) ~ (psp\0\p sp ) + - / dpdh p\ 1 (p)p sp {h) 

th 1 J -oo 

x (p sp \0\p sp , {h,p})e~ itu ^ Pap ’^ h ' p ^~ Ss ^ Pap ^ h ’ p ^ + mirr. 

( 11 ) 

Since the dispersion relation, as well as the differential 
overlap coefficient Ss[p, {h,p}], splits in terms of parti¬ 
cles and holes u>\p, {h,p}] = u;[p,p] — w[p, h], the inte¬ 
grals can be approximated by evaluating each of them 
in the saddle point of the single-particle dispersion rela¬ 
tion w[p, A] which for any smooth distribution p(A) is in 
A = 0. Therefore if k is the order of the first non-zero 
derivative in A = 0 of e a 0 A) ; the approach to the steady 
state value of all the local operators with a finite expec¬ 
tation value on the saddle point state ((p S p\0\Psp) ^ 0) 
is given by a power law as follows 

Aa(tW (fe+2) Vfc odd , A^t) ~ t~ ( ' k+1 ' 1 Vfc even , 

(12) 

where A ^(t) = ( 0{t)) — lim t _ > , oo (0(t)). The power law 
decay of correlations is a consequence of the creation, by 
the quench, of a finite density of holes around A = 0, 
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FIG. 1. (a): (color online) Time evolution of 52(2 = 0, f) as a 
function of time for different values of the post-quench interac¬ 
tion 7 = 4, 8 , 16, (from top to bottom) in the thermodynamic 
limit with fixed density n = 1. The data are obtained using 
equation |9]l and by averaging over 10 finite size realizations of 
the saddle-point states with a system size of N = 96 particles. 
The lines on the right respectively indicate the steady state 
values in the thermodynamic limit as given in [SB], (b): (color 
online) Time evolution of gi{x = 0) as a function of time for 
7 = 4 and different system sizes: N = 6,8 (red, blue line) 
and in the thermodynamic limit (black line). The finite size 
asymptotic values (red and blue lines on the right) are shown. 
The data for N = 6 and N = 8 are obtained by performing 
the full double sum over the Hilbert space (|2| and dividing by 
their initial value 52 ( 2 : = 0,t = 0) = (1 — jj). Their asymp¬ 
totic values correspond to the value of the diagonal ensemble 

E X e (A|:(p(0)/n)2|A) 

g 2 (x=0,t=0) 


g 2 (x = 0 ,£) 



giving a finite density of states for small-energy (zero ve¬ 
locity) particle-hole excitations in this region [50] (See 
figure [ 2 ] panel (b)). Therefore the contribution of the 
power law is proportional to the density of holes around 
A ~ 0 in the post-quench saddle point state which is 
large for distributions with a large (extensive) entropy. 
Any initial state with an extensive amount of energy 
(’I r o|.H'| 1 &o) ~ Leo shows therefore a relaxation as power 
law although its contribution to the whole time evolution 
becomes less and less visible as eo decreases. Note that 
up to now we assumed that the operator O conserves the 
total number of particles. For operators adding (or re¬ 
moving) one extra particle to the system the power law 
is simply replaced by f - ( fc + 2 )/ 2 (t~( k + 1 )/ 2 for k even). 

Finally it is important to note that the same large time 
decay is expected also for systems with bound states (un¬ 
der the string hypothesis) [39MT] . However in this case 
the full time evolution from t = 0 + can only be recov¬ 
ered by including other classes of high energy excitations, 
namely recombinations between bound states of different 
masses (strings of different lengths [BUI). 

Time evolution in the interacting Bose gas As a spe¬ 
cific example of the general method, we now focus on the 
Lieb-Liniger model for a S —interacting Bose gas, defined 
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FIG. 2. (a): (color online) Log-log plot of the time evo¬ 

lution of A 92 (£) for different values of the post quench in¬ 
teraction 7 = 4,8,16 (from top to bottom). The red lines 
f(t) = const( 7 ) — 31ogt are guide for the eyes showing the 
approach to the equilibrium value as ~ t~ 3 for all the con¬ 
sidered values of the post-quench interactions, (b): (color 
online) Schematic representations of the most relevant parity 
invariant particle (filled dots) with its respective hole (empty 
dot) excitations on the saddle point filling function (here for 
7 = 4). For small times the high energy excitations corre¬ 
sponding to particle-holes along the tails (blue dots) need to 
be included. For large times on the other hand the relevant 
excitations are particle-holes deep in the center of the distri¬ 
bution A ~ 0 (red dots) which are the ones responsible for 
the power law decay of the correlations towards their steady 
state values. 


by the Hamiltonian [Mj (setting h = 2m = 1) 


j>k 

The initial state is chosen to be the ground state in the 
absence of interactions 7 o = 0 , where 7 = c/n effec¬ 
tively parametrizes the coupling in the thermodynamic 
limit. This state is known as the Bose-Einstein conden¬ 
sate (BEC) state |BEC) and it is spatially structureless 
in all coordinates, (x|BEC) = >y/ 2 - We consider the 
post-quench time evolution of the static density moment 
g 2 {x = 0), measuring the rate of two-body inelastic pro¬ 
cesses in the gas [65] which can be experimentally ac¬ 
cessed through the measurement of the photoassociation 
rate [57] 

g 2 (x = 0 ,t) = (BEC\e zHLLt :(/5(0)/?r) 2 : e~ iHLLt \BEC) . 

(14) 

The density operator is defined as p( x) = (a:)(a:), 

where the bosonic operators ’if'(a;), iIU(:r) satisfy the 
canonical commutation relations [\t(a;), ^(a/)] = S(x — 
x'). The overlaps and in particular the generalized one- 
particle overlap coefficient have been computed in [38] 

So EC (A)=log^^ + ^ , (15) 


J^S(xj-Xk). (13) 


N d 2 

H ^ = -T,Yx f + 2c 

3 = 1 3 


where the branch-cut of the logarithm is chosen such that 
So EC (—A) = —Sq EC (A). The one-particle energy and mo¬ 
mentum are given by fco(A) = A and eo(A) = A 2 . The 
function s EEC (A) determines the saddle point state which 
can be analytically written in terms of Bessel functions 
of the first kind I n (z ) [381 


Psp (A) 


a sp (A) 


7 1 da sp (A) 

47 rl + a sp (A) d 7 


The matrix elements between the eigenstates of the 
model are given in mm- The sum is performed by av¬ 
eraging over different finite size realizations | A sp ) -7 \p sp ) 
of the saddle point state and evaluating the relevant ex¬ 
citations via an adaptation of the ABACUS algorithm 
[68U70] to generic highly-excited states. In figure [I] the 
time evolution computed via the quench action approach 
shows that even for values of the coupling constant 
that are far from the two perturbative regimes (weak and 
strong coupling) we recover the initial BEC value of the 
correlation lim t h<72(^ = 0 ,t = 0 + ) = <72(2 = 0 )bec = 1 
[65l . The thermodynamic results allow to extract their 
large time decay to their steady state values as in figure 
[ 2 ] This follows the expected t~ 3 law which is a conse¬ 
quence of (12 1 and of the behavior of s EEC (A) around 
A = 0 


exp (s EEC (A = 0)) = 0 


dexp (s EEC (A)) 


dX 


A —0 


7^0 


V 7 > 0, 
V 7 > 0 . 


(17) 

(18) 


This shows that the relaxation following a power law is 
present for any post-quench coupling constant 7 , even 
in the limit of small interactions. This is in contrast to 
the predictions of the Bogoliubov approximation where 
the decay is predicted to be exponential for small 7 m 
Note that the behavior of the overlap as in is also in¬ 
dependent of the initial value of the coupling constant. It 
is related to the fact that for quenches from the ground 
state of the theory with a coupling 70 > 0 to the gas 
with a finite coupling 7 > 0 the eigenstate with the max- 
irnal overlap e~^>- is clearly the ground state of the final 
theory. This leads to the divergent behavior of the gener¬ 
alized single-particle overlap for small values of the rapid¬ 
ity, limA-j-o e _s »°A) 00 which leads to Therefore 

the same power law is expected for any interaction 
quench y 0 —> 7 > 0 inside the repulsive regime of the one- 
dimensional Bose gas (for quenches to the free bosonic 
theory 70 —1 7 = 0 see mm- 

Conclusions We showed how the quench action ap¬ 
proach allows to reconstruct the whole post-quench time 
evolution of an integrable system from data contained in 
the thermodynamically leading part of the overlaps. In 
particular we presented an argument to predict the power 
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law behavior for the late times approach to equilibrium of 
local observables. This is a direct consequence of the cre¬ 
ation in the gas of macroscopic excitations with vanishing 
velocity which is a generic feature of the model itself, in¬ 
dependently of the quench protocol. The question if an 
adaptation of the non-linear Luttinger liquid approach 
for equilibrium correlations j48H50| can be implemented 
to compute the late time dynamics after a quench will be 
addressed in forthcoming works. 

As a proof of principle we computed the time evolution 
in the Lieb-Liniger model of the static density moment 
g 2 (x = 0,t) after a quench from the Bose-Einstein con¬ 
densate. This represents a rare example of a full post- 
quench time evolution of a truly interacting model and 
therefore it can be directly connected to experimental 
results in ring-likc geometries EU, box-like potentials 
m or any other experimental realization of the one- 
dimensional Bose gas where the confining trap influences 
time scales which are much larger than the relaxation 
time of one-point functions as g 2 (x = 0,f) [S', ,T2]. The 
comparison between the finite size calculations and the 
thermodynamic limit in figure [l] shows indeed that for 
short times the relaxation processes are well approxi¬ 
mated by N ~ 10 particles. This also underlines the 
importance of obtaining exact results in the thermody¬ 
namic limit, that can be used to test numerical simula¬ 
tions for small system sizes as done in [SB]. The method 
can be extended to two-point functions as the dynami¬ 
cal density-density correlations of the gas and to other 
models as the XXZ spin chain [39, l4Tj . 
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